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We studied the process of decoherence induced by the presence of an environment in acoustic black 
holes, using the open quantum system approach, thus extending previous work. We focused on the 
ion trap model but the formalism is general to any experimental implementation. We computed the 
decoherence time for that setup. We found that a quantum to classical transition occurs during the 
measurement and we proposed improved parameters to avoid such a feature. We provide analytic 
estimations for both zero and finite temperature. We also studied the entanglement between the 
Hawking-pair phonons for an acoustic black hole while in contact with a reservoir, through the 
quantum correlations, showing that they remain strongly correlated for small enough times and 
temperatures. We used the stochastic formalism and the method of characteristic to solve the field 
wave equation. 
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I. INTRODUCTION 

Hawking effect, i.e. the particle creation process that 
gives rise to a thermal spectrum of radiation outgoing 
from a black hole [1] , is a prediction of quantum field the- 
ory in curved space-time. This effect together with its en- 
tropy complete the interpretation of black holes as ther- 
mal objects. On one hand, it is believed that the heart of 
a theory that unifies quantum mechanics and gravity lies 
in understanding the nature of this thcrmality. On the 
other hand, it is also important to collect experimental 
evidence in order to gain insight into this phenomenon, 
but this is practically impossible since black holes' tem- 
peratures are less than nK. Since the phenomenon has 
not been observed experimentally, it is of crucial impor- 
tance that all the assumptions underlying the Hawking 
effect be carefully studied so as to try to understand the 
process. 

W. G. Unruh proposed an analogue gravity hydrody- 
namical model where phonons propagate in a fluid with 
a subsonic and supersonic regime [2]. This model obeys 
the dynamics of a massless scalar field near a black hole 
and provides a possible experimental implementation to 
study the Hawking effect. Subsequently, there have been 
several more realistic proposals that involved BEC [3], 
moving dielectrics [4], waveguides [5], slow light systems 
[6], among others. 

On general grounds, an acoustic black hole is a system 
in which phonons are unable to escape from a fluid that 
is flowing more quickly than the local speed of sound. 
Therefore, trapped phonons are analogous to light in real 
gravitational black holes. The phonons propagating in 
perfect fluids exhibit the same properties of motion as 
fields, such as gravity, in space and time. 

The current proposals do not provide conclusive evi- 
dence of the Hawking effect (see for example [7]). Nev- 
ertheless, we believe that a particular one by Horstmann 
et al., [8], provides a promising setup. This system con- 



sists of a circular ring of trapped ions moving with an 
inhomogeneous velocity profile emulating a black hole, 
as explained above. The signature of this quantum radi- 
ation is the correlation between entangled phonons near 
the horizon, [9], which can be measured by coupling the 
ions' motional degrees of freedom to their internal state, 
[10]. 

Given that we are interested in the quantum nature of 
the effect, in this paper we make a study of the interac- 
tions of the acoustic black hole with an environment and 
how this induces decoherence in the system. 

The presence of an environment can destroy all the 
traces of the quantumness of a system. All real world 
quantum systems interact with their surrounding envi- 
ronment to a greater or lesser extent. As the quantum 
system is in interaction with an environment defined as 
any degrees of freedom coupled to the system which can 
entangle its states, a degradation of pure states into mix- 
tures takes place. No matter how weak the coupling 
that prevents the system from being isolated, the evo- 
lution of an open quantum system is eventually plagued 
by nonunitary features like decoherence and dissipation. 
Decoherence, in particular, is a quantum effect whereby 
the system loses its ability to exhibit coherent behavior. 
Nowadays, decoherence stands as a serious obstacle in 
quantum information processing. All in all, it should be 
crucial to control decoherence in order to plan a concrete 
setup for measuring the Hawking radiation as a quantum 
feature of analogue gravity models. 

In this paper, we show a more detailed and deeper 
treatment of the problem given in [11], and we present 
new results concerning the evaluation of the decoher- 
ence time and analytical expressions for correlation func- 
tions of the quantum fields. We continue using the field- 
theoretical description in order to present a derivation 
applicable to any implementation of an acoustic black 
hole. Nevertheless, given our interest in the ion trap we 
provide numerical results corresponding to that particu- 
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lar setup. 

The article is organized as follows: in Sec. II we show 
the specific model we are using. In Sec. Ill we present 
the model of environment and show the non-equilibrium 
dynamics of the system. Sec. IV contains the estimation 
of the decoherence time and we show, in Sec. V the 
dynamics of the entanglement through the analysis of 
correlations. Finally, we summarize our results in Sec. 
VI. 



II. A PARTICULAR ACOUSTIC BLACK HOLE: 
THE ION RING 

As stated above, we will work in the context of the field 
theoretical description. Given that we are interested in 
the ion ring model, and with the purpose of introducing 
the subject of analogue gravity with a specific example, 
we will explain in detail this set up from first principles. 
We will begin with the discrete description of N ions of 
mass m in a circular trap of radius R to end with the 
description of a massless scalar field in a curved back- 
ground. 

Following the study done in Ref. [8], the Hamiltonian 
describing the ions of the circular trap is given by 



where we call r = 0.05 • T the collapse time, which is 
much smaller than T, and after that time the asymptotic 
profile is given by the parameter v m ; n (i ^> r) = w m ; n . It 
is important to notice that the measurement can only be 
performed during one period T, since after that time the 
system's classical configuration loses its stability, see Ref. 
[8] . Therefore, we assume that any possible measurement 
ought to end after one revolution, but must last longer 
than r in order for the acoustic black hole to form, as we 
will see below. 

The problem can be linearized for small perturbations 
of the trajectories 9t = 9® + 89i as 
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where V e {0i,i) is an external field potential correspond- 
ing to an electric field that induces the classical trajec- 
tories 9®{t) and V c (9i, ...,9jy) is the Coulomb potential 
between the ions. Those are such that the velocity profile 
is 



v(9,t) 



a 



71 



= < 



^max 
P-l 
^min 



j-27T+g£ 
72 



< 9 < 9 H -7i 

-7i < (0 - H ) < 7i 

9h + 7i < 9 < 2tt - 9 H - 72 

-72 < (9-2t: + 9 h ) < l2 

2tt - 9 H + 72 < 9 < 2tt, 



where /3 = (w max + v min )/2 and a = (v max -v min )/2. The 
minimum and maximum velocity are constrained by the 
condition that each ion has to make one revolution dur- 
ing a period T. It is important to notice that we use an 
approximate velocity profile, as introduced in Ref. [11]. 
The real, as explained in Ref. [8], must be C 3 but ours 
is a good approximation appropriate to our calculations. 
Taking this into account, the parameters 7, and 9r do 
not exactly match those of Horstmann et al. The system 
may be prepared in an initial thermal state with temper- 
ature To. An illustration of the setup and the velocity 
profile can be found in Figs. 1 and 2 of Ref. [8]. 

Initially, the velocity profile is constant with v m - m (t — 
0) = v max (t = 0) = 27r/T but during a time r the profile 
change such that u m i n (i r) < 2n/T and v max (t » 
r) > 2tt /T in the following way 
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We are interested in the continuous description of this 
system, described by the field $(6» = 0?(t),t) = 59i{t). 
After diagonalizing the previous hamiltonian and going 
to the continuous limit we end up with the Lagrangian 
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where D(9, k) = c(9)k + 0(k 3 ) is the dispersion relation. 
The speed of sound is given by 
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where n is the local density of the ions and Q their elec- 
tric charge. The conformal factor is p{9) — mR 2 n(9) = 
mR 2 (N/(v(9)T)). 

In Ref. [12], Unruh has shown that the Hawking ra- 
diation is robust against small deviations from a linear 
dispersion relation, so we take simply D(9,k) — c(9)k. 
Moreover, the conformal factor can be included to first 
order in the definition of the field, redefined as = 
Taking this into account, we can describe the field by 
means of an action written in the suggestive form (this 
description is the starting point also in Ref. [11]) 
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where x M = (i, x = 9) and the effective metric is 

ds 2 = (c 2 - v 2 )dt 2 + 2vdxdt - dx 2 . (7) 

Of course this system is non-relativistic and it must be 
thought of as an effective sigma-like model describing the 
phonons in the ion trap. Nevertheless, this is indeed an 
analogue black hole, since it presents an event horizon in 



the points where v = 
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This analogy can be clarified 
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then the metric is 

ds 2 = (c 2 - v 2 )dr 2 



dx 2 
1 — v 2 /c 2 ' 



(8) 



and expanding linearly around the angle 9h of the event 
horizon where v(6u) = c(9h), we recover the Schwarchild 
metric. This field 4>{x) describing the phononic exitation 
of the continuous array of ions can now be quantized. Re- 
peating the Hawking radiation derivation, the Hawking 
temperature can be obtained 
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As suggested by Balbinot et al. [9] and studied by 
Schiitzhold and Unruh [13], the signature of the quantum 
radiation, in contradistinction with a stimulated emis- 
sion, and what would be measured in the proposal, is 
the peak present in the correlation (5pl(0)Spl{0')) , be- 
tween two points inside and outside the acoustic black 
hole. The magnitude 5pL is the left moving component 
of the canonical momentum conjugate to SO. 

This feature reflects the entanglement between the 
Hawking pair of phonons emitted near the event hori- 
zon. Following the discussion in Ref. [9], this can be 
seen writing the "in" vacuum in terms of the "out" vac- 
uum as a squeezed state 
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where a( esc )>v tr )t are creation operators for, respectively, 
the outgoing escaping and trapped modes. 

This magnitude was calculated numerically and the re- 
sults are shown in Fig. 9 of Ref. [8] . When one increases 
the initial temperature Tq of the system, the entangle- 
ment is lost due to a quantum to classical transition in- 
duced by thermal effects. In [8] they computed, as a mea- 
sure of the entanglement, the logarithmic negativity, see 
Fig. 14 of Ref. [8]. On one hand, for fixed temperature, 
the system needs a certain amount of time to generate 
the entanglement through the creation of the Hawking 
pairs. As the temperature increase above a threshold, 
~ 100 -Tjj, the entanglement is lost and the system starts 
behaving classically. This is also reflected in the corre- 
lation at the same temperature scale. In this paper we 
will be interested in the quantum to classical transition 
induced by the non-equilibrium dynamics with an envi- 
ronment instead. 

Another feature that can affect the measurement is 
the presence of a noise in the force used to produce 
the trajectories. This stochastic component of the force 
can be described by the parameter 7 defined such that 
(noise) = 7 x (mean force). The requirement that the 
peak in the correlation is distinguishable imposes the fol- 
lowing bound 



Although the problem can be treated in a discrete fash- 
ion, we choose to stick to the field description since the 
action given in Eq. (6) is common to every realization of 
acoustic black holes, not restricted to ion traps. For ex- 
ample, in the case of a BEC acoustic black hole it emerges 
from the Gross-Pitaevski equation or in the waveguide 
setup from Maxwell's equations, etc, but the point is that 
the system must always be described with an action as 
Eq. (6) in some regime. Therefore, it is important to 
notice that the analysis carried in the following sections 
can be applied to any set-up. 



III. CHARACTERIZATION OF BOSONIC 
ENVIRONMENTS AND THE 
NON-EQUILIBRIUM DYNAMICS 

A. Interaction Model 

Our aim in this paper is to study the behavior out of 
equilibrium of this system while interacting with a quan- 
tum environment, since this coupling acts as a mecha- 
nism that induces decoherence and the system starts be- 
having classically after the decoherence time. As usual 
for this kind of tasks we use the Schwinger-Keldysh or 
"in-in" formalism. 

As explained in the previous section, our system is de- 
scribed with the action of a massless scalar field in a 
dynamical background, 



So[<t>] = \ 



(12) 



Following the quantum Brownian motion (QBM) 
paradigm, the environment is described by a continu- 
ous array of bosonic quantum harmonic oscillators dis- 
tributed in each position of the circular trap, following 
Ref. [11]. The bath is at rest with respect to the labora- 
tory and it will be represented by the following degrees 
of freedom q v (6,t) with the action 
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dvl{y) I d 2 x[ql{x,t)-v 2 ql{x,t)]. (13) 



The function I{v) corresponds to the mass of each os- 
cillator that compose the environment. The interaction 
between the system and the environment is given by the 
following term in the action 

Sint[(f>,Qv] = -J dv J C(u)d 2 x 4>(x)q u (x), (14) 



in such a way that the total action is 

S[(j),q v ] = S [4>] + S £ [q, y } + S in t[^,q v ] 



(15) 
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The nature of the environment depends strongly on 
the specific model of black hole considered. For example, 
in the ion ring proposal the velocity profile is generated 



4 



by the action of an electric field, which in turn is pro- 
duced by plane-parallel electrodes. The irregularities of 
its surface produce fluctuations in the force they induce. 
The noiseless component of this force is included in the 
effective action of the system in Eq. (12), such that the 
analogue gravity model is achieved. The pure noise com- 
ponent due to the fluctuations can be introduced as an 
environment and the oscillators represent the correspond- 
ing modes of the stochastic electric field. The details of 
the possible natures of an environment was extensively 
discussed in Ref. [14]. The fact that we are thinking of 
an electric field coupled to the ions' coordinates justifies 
the bilinear coupling used here. Nevertheless, in other 
analogue gravity models as long as the environment is 
bosonic, this model is fairly general. For example, if the 
field is coupled to the ions' velocity instead of their posi- 
tion, then this amounts only to a change I(v) H ► I(y)v 2 , 
as can be seen below. Therefore, the information of each 
particular case is encoded in the spectral density. 

To simplify the analysis we take an ohmic environment, 
although more general environments do not change sub- 
stantially the results as can be seen in Ref. [15] in the 
context of quantum Brownian motion. 

In Ref. [16] it was shown that the presence of an en- 
vironment with scales bigger than the Planck scale can 
generate instabilities such as Miles-type instabilities that 
jeopardize the detection of the Hawking radiation as a 
quantum induced effect. The environment is at rest with 
respect to de laboratory and within the approximation 
used, both the environment as the ions are in the linear 
dispersion relation regime, far from being in the Planck- 
ian regime, where the correct description is the Coulomb 
chain expression. Therefore, the model used here does 
not present any instability whatsoever induced by the 
environment. 



In the following we will present the formalism used 
to study the non-equilibrium dynamics of this composed 
system. 

B. Reduced Density Matrix 

The dynamics of any system out of equilibrium is de- 
scribed by its density matrix together with its evolution 
in time. The matrix elements of this operator are given 
by 

p{ct>,q;<t>',q'\t) = (ct>,q\p{t)W,q')- (16) 

We use an uncorrelated initial state between the system 
and environment, in such a way that 

p{U)= ps(ti)®pe(ti). (17) 

These initial density matrices correspond to a thermal 
state with temperature To = (/cs/3) _1 , i.e. p ~ e~^ H . In 
principle the initial temperature of the system and the 
environment can be different and the result we obtain 
below for the decoherence depends only on the initial en- 
vironmental temperature. Nevertheless, it is more appro- 
priate in the experiment to think of an initial equilibrium 
state between the system and the bath. 

More general states do not change substantially the 
process of decoherence, as can be seen in Ref. [17]. 

The reduced density matrix, which represents the con- 
cept of coarse graining the environmental degrees of free- 
dom, which are of no interest to us, is defined in the usual 
way as the partial trace over their degrees of freedom 

p r (<t>,<fi'\t) = Tr e p. (18) 



The generalization of the non-equilibrium dynamics of a quantum harmonic oscillator to the study on fields is 
straightforward and can be found in Ref. [18]. The evolution in time of the reduced density matrix is given by the 
following expression 

p r (<f>,<f)';t) = J VfcVtfiJrty, <j>'\t\<l>i,(l>'tiU)ps{fc,<t>'i\k), (19) 
where we defined the evolution operator associated with the reduced density matrix 

JrM-AfcAlU) = [ 2ty r 2y e i < s °M-*°l*'lV\F[^']. (20) 
Jtpi Jtf>'. 

In this expression the influence functional is defined following the Feynman- Vernon treatment 

T[4> ) 4f] = ^ anf{ *'*' A,n = I Vq v ps(qi,4)exp!i(S e [q v }-Se[q' v } + S int [^q]-S int [^,q])/n\, (21) 

JCTP I J 

where the prime and unprimed represent both branches of the closed time path curve over which we integrate, see 
Ref. [18, 19]. 

Therefore, the evolution of the system is dictated by the coarse grained effective action 

S cS [0, 0'] = S [<t>] - So [<j>'\ + Sif [0, 4>'\ • (22) 
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The computation of the Feynman- Vernon influence action is identical to the QBM case and the result can be cast 
into the form 



S 



IFl 



} = / d 2 xd 2 x'(b-(x) [B(x,x')^ + (x') + ^N( 



V) > 



(23) 



with the usual definitions 



b')/2. The kernels D(x, x') and N(x,x') have the same expression 



that the QBM case regarding the temporal behavior, and they are local in space. 

/>oo 

D(x,x') = / duJ(u)sinv(t-t')Q(t-t')S(x - x'), 
Jo 



N(x,x') = 

where the effective spectral density J{v) is defined as 



dvJ(v) coth — cos v{t — t')S(x — x'), 



J(iz) 



As we stated previously, we consider an ohmic environ- 
ment and this translate into the following form of the 
spectral density 



(27) 



where 7 plays the role of an effective coupling constant 
and f{y) is a generic cut-off function whose effect is to 
regularize the expression, and to this aim it must satisfy 
the following requirements 



f(v = 0) = 1 and j(y > A) 



0. 



(28) 



Starting with the effective action S c s, we can write 
the semiclassical equation of motion, analogous to the 
stochastic Langevin equation for the field, i.e. 



1 d_ 



V99' 



dx v 



+ / dsT>(t,s)(j)(s,x) = £(x,t). 



The field £ is a stochastic force with a gaussian proba- 
bility distribution with zero mean (£) = and two-point 
function given by 



{Z(x)Z(x , ))=hN(x,x'). 



(29) 



This stochastic force can be identified as the noise pre- 
sented in the previous section, quantified with the pa- 
rameter 7. 



C. Estimation of the Effective Coupling Constant 

To make contact with the parameters used in Ref. [8] 
we will calculate the relationship between the effective 
coupling constant 7 coming from the environmental spec- 
tral density and the parameter that characterizes the 
noise, the fluctuations in the force applied to the ions, 
introduced above as noise = "fF. 



i( v y 
I 



(24) 
(25) 

(26) 



Making a comparison between the equation that de- 
fines 7 with the Eq. (29) it is possible to identify the 
stochastic force with the fluctuations jF. Using the 
equation D<f> + ...=£ one can obtain mU[R50) + . . . = 
mi?£ / y/p and therefore 



mR . 

7^ 



£rms = lF, 



(30) 



where £ rms = yj is the root mean square of the 
stochastic force. 

The mean force can be estimated as F ~ rnR(v mayi — 
Vmin)/ T - Therefore, coming back to the previous expres- 
sion 



4=^ = 7"^ 



(31) 



The noise kernel is given by = fiN(O) and if one 
uses a cut-off function of the form f(u) = e~ v l K we ob- 
tain 



N 



Therefore £ r 



1 



1 



7 = 7" 
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/ H/2 ; yA, and 
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(32) 



(33) 



If we assume that r is the minimum amount of time 
during which the distribution changes appreciably (as we 
did to compute the mean value of the force) and it is 
of the same order that the magnitude of the time scale 
associated to the environment, then the Lorentzian cut- 
off frequency satisfy the relation At rj 1 and we finally 
arrive to the desired expression 



7 = l—^{ v n 



(34) 
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In accordance to the results presented in the Ref. [8] 
the bound of the coupling constant proposed originally is 
7(7) < 5- 10~ 6 . As explained in the previous section, this 
bound comes from the requirement that the fluctuations 
in the trajectories due to the noise in the force do not 
wash out the characteristic peak in the correlation. 

In the following section we will learn that this bound 
is not appropriate since in this range decoherence occurs 
even before the acoustic black hole is formed, at a time of 
the order of r. To increase the decoherence time to values 
bigger than the measurement time, of order T, the mag- 
nitude of the noise present in the system, characterized 
by 7, must be appropriately reduced. 



IV. DECOHERENCE TIME 

A. Master equation for the reduced density 
matrix. 

To obtain the decoherence time first one has to obtain 
the relevant coefficients included in the master equation. 
In particular, those corresponding to diffusive effects. In 
order to obtain this equation, the procedure consists of 
deriving the evolution operator with respect to time to 
generate the derivative of the reduced density matrix. 
The next step is to multiply 3 r p r {ti) and finally integrate 
the initial conditions thus generating p r (t) instead of its 
initial value [18, 19]. 



The derivative of J r is 

d l ( d 1 

ifr-Qj.Jr[4>f,<l>f,t\4 > i,<l>i,ti\ = Jr[4>f,(l>f,t\4>i,<t>i,0]s ~ -K^S e s[(t>ch4 > ch t ] h ( 35 ) 

and using the following identity that comes from the causality of the dissipation kernel 

Sw[<f),(f)',t + dt}= SiF[<f>,(j)',t]+dt J dx4>~(x,t) J ds{T)(t,s)(t) + (x,s) +iN(t,s)</>~(x,s)}, (36) 
one can obtain, after performing the functional integrals 

^ Pr (frc( > > l t) = ~(<t>\[Hs,Pr(tW}-l J dx{<t>(x)-<t>'{x)) j dt' 

X |n(*, [$ - *'] (0, 0'; *') - l -B(t, t') [* + &] (0, <f>',t)\. (37) 

The first term that comes from deriving the free part of the effective action generates the usual Liouville-von Neumann 
term. The second contribution was defined as 

0= / VxVx' e iSrf <V* Pr (^ O; ^ ; o)0(x,O, (38) 

J <p(t,x)=4>(x),4>'(t,x)=(j) l {x) 

and a similar expression for <&'[<?!>, <f>'; t'\ with <j)(t') i-> 4>'(t'). In the case of weakly coupled environments, 7 ~ 10~ 6 
and since this terms are already beyond linear order with respect to the coupling constant because of the presence 
of the kernels, the expression above can be approximated by S e g[4>, <f>'] ~ Sq[4>] ~ Sq[4>'] in the functional integral. 
This way, we can see that this functions are the matrix elements of the operators $(x,t') = 4>(x,t' — t)p r (t) y 
<&'(x, if) — p r (t)4>(x,t' — t). Using this expression we can finally obtain the master equation from 

H^- t p r (t) = -i[H s ,p r (t)}-J^ dTdx{N(t,t-T)[$s(x),[fa,-T),fir(t)]]^^ (39) 

where 4>s{x) is the field operator in the Schrddinger picture and <J){x,—t) corresponds to the Heisenberg one. To 
study the decoherence, we are only interested in the term proportional to the noise kernel since those are the ones 
that generates the necessary diffusion terms in the master equations. To achieve this, we need to replace the solution 
of the Heisenberg equations of motion (EOM). 



To continue we will obtain the solutions of the Heisen- be cast in the following form 
berg EOM, which coincide with the classical EOM of the 
field. Since we are in 1 + 1 dimensions, the solution can 



= /(u)+ 5 (v), 



(40) 
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where we define the null coordinates associated with the 
effective metric 

dx 



u=t- 
v = t + 



c{x) + v{x) 
dx 

c(x) — v(x) 



Xu 

= t — Xy. 



Therefore, we can study the decoherence for each mode 
u and each mode v separately 



(x,t) = u / v coscj(i- x u/v ). 



(41) 



The allowed frequencies u> can be found requiring the fol- 
lowing condition since the spatial dimension is compact, 
x = x + 2n7r, and one has to impose periodic bound- 
ary conditions over the ring cosw(t — x u / v (x + 2nir)) = 
cosw(i — x u / v (x)). We take as the maximum allowed fre- 
quency to be w max ~ N/T. The solution is continuous 
in the subsonic-supersonic transition regions. The con- 
tribution proportional to the canonically conjugate mo- 
mentum II5 can be discarded since it does not generate 
a diffusive term. 

We will study separately the decoherence present in 
each mode of frequency u) and both u and v modes. In- 
serting this solution of the EOM in Eq. (39) the relevant 
term is given by 

drdxN(t, t — r)[(f) u / v cosu>ar u / v , 

[4> u/v CO8 0j(T + X u / v ),p r (t)]\. (42) 

Taking the appropriate expectation value one obtains 



(0u/ v -0u/ v ) 2 / drdxN(t,t - t)cosux u/v 
Jo 

x coso;(a; u / v + r)p r {<j>, <j>') 

(0 u /v - 0u/ v ) 2 / drdxN(t, t — r) cos uix u / v 
Jo 

)swi u /, coswr — s'mu!x u / v sinwr^j p r ((j>, 4>). 



To estimate the decoherence time we study trajectories 
with (4> u /v ~ ^u/v) 2 ~ pfi* ! ) where 5 = 2ir/N, the mean 
separation of the ions since it is the smallest length scale 
of the system. We use the usual definition of the normal 
and anomalous diffusion coefficients respectively 

D(t)= [ drN{t,t-T) cos lot 
Jo 

f(t) = I dTN{t,t-r)smuJT, 
Jo 

and we also define the following coefficients 

r 2ir 

^2 



(43) 



Vi u / V 



^2u/v = 







dx(cosu)x u / v ) 



2tt 



dx cos ujx u / v sin ujx u / v . 



(44) 



Therefore, the diffusion term we obtain is given by 

h Pr ~ ~ P S 2 {D(t)V lu/v + f(t)V 2u/v } Pr . (45) 

This term contributes in the following way to the solution 
of the master equation 



Pr 



p^exp 



\p I di><D(l')V lu v - f[l')V 2u v 



(46) 

where p u represents the unitary evolution of the reduced 
density matrix. The decoherence time is defined based 
on the following relation 



P 



dt'(D(t , )V lu/v + f(t')V 2u/v )8 2 



1. 



(47) 



This corresponds to the fact that for times t > tjj the non 
diagonal elements of the density matrix that encompass 
the quantum coherence effects are suppressed. In order 
to compute it, we will need explicit expressions for D(t) 

y fit). 

Lets start with the normal diffusion coefficients at zero 
temperature Tq = 



D(t) 



7 2 



dv I ds — 
11 Jo 1 



2 cosi/scosws, (48) 



where the expression was written with a Lorentzian cut- 
off to attenuate high frequencies. According to Ref. [20], 
the explicit calculation of this integral gives the following 
result 



D(t) = 



I 2 



Shi(At)l — cos ujt cosh At 



- sin jjl si nli \l j — (3hi(At)( — cos u>t sinh At 



- sin u)t cosh At 



Si(wt) 



where Shi(a:) and Chi(x) are the hyperbolic sine and co- 
sine integral respectively and Si(x) is the sine integral. 

On one hand, for times much larger than the frequency 
cut-off scale, i.e. t 3> A -1 , this expression can be approx- 
imated as 



D(t) 



I 2 



if) 



rSi(wt). 



(49) 



On the other hand, for times much larger than the scale 
imposed by the field frequency, t 3> the sine integral 
can be approximated as Si(wt) « 7r/2 and for frequencies 
much lower than the cut-off uj < A, 1/1 + (^) sa 1. 
Finally, one obtains the expression 

D(t)*f^. (50) 

Therefore, the integral of the coefficient is given by 



dsD(s) 



7 1 ut D- 



(51) 



The same analysis can be done for the anomalous dif- 
fusion coefficient and can be found in [20] . The result for 

t > w -1 , A -1 is 



/(*) 



7 2 , A 



(52) 



As in the previous section, we can estimate the charac- 
teristic time of the environment as A -1 ~ r, then after 
performing the integral, the anomalous diffusion coeffi- 
cient can be estimated as 



7 2 1 
dsf{s) « — Trlog — uit D . (53) 

Z LOT 



Putting the previous results together, for an initial state 



0.006 




(a) 




(b) 



FIG. 1. Plot of the functions (a) log (ojt)~ 2 V 2u / v < V lu / V , 
showing that the anomalous component is negligible; and (b) 
Vi u /V2v with respect to the mode frequency to, showing that 
the decoherence time is the same for both u and v modes. 



with zero temperature Tq — 0, the decoherence time is 



given by the following expression 
t D (T = 0) 



4/j 



P S 2 j 2 ttlo(V 1u/v + 2 log ^V 2u/v ) 
_ 2h 2 
r y 2 AvS 2 uJTrp 2 

x(V lu/v + 2\og^-V 2u/v )-\ (54) 

where Av = v max — i> m in- If we focus in the ions' ring and 
the experimental parameters proposed for the system, 
then the expression can be simplified even more, since 
we are in a possition to approximate 



2 log 



I 



^2u/v<^iu/v and P lu ~P 2v = P, (55) 



as can be verified in the plots shown in Fig. 1. Finally, 
we obtain the following result for the decoherence time 



t D (T o = 0) 



2h 2 



(56) 



7 2 Av<5 2 u;7rp 2 V 
In Fig. 2 the result for the decoherence time as a function 




l.xio- 6 2.x 10 



3.X10- 6 4. xlO- 6 5.x 10" 



FIG. 2. [Color Online] Dependence of the decoherence time 
as a function of the relative noise in the force 7. The colored 
band includes all the allowed frequency range of the solutions. 

of the force relative noise 7 is shown. The decoherence 
time is much shorter that T for the bound presented in 
the Ref. [8], 7 ~ 5 • I0~ 6 . As explained in section II, 
the period T is also of the same order of magnitude as 
the time needed to perform the measurement. For times 
larger, the system becomes classically unstable and for 
times smaller the black hole would have no time to de- 
velop. The system, for the experimental parameters pro- 
posed in Ref. [8], shows decoherence in a time-scale too 
short and this would make impossible the measurement 
of the aspects of the Hawking effect, which is purely of 
a quantum nature. Moreover, decoherence happens in 
a time of the same order of magnitude as the collapse 
Id ~ T i therefore there is no time for the acoustic black 
hole to form. 
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To find a solution for this issue, we impose a decoher- 
ence time appropriate to the experiment, to be specific, 
no shorter than the following bound > 100 • T. 

For this to be satisfied and to avoid the quantum 
to classical transition induced by the environment, the 
bound for the coupling with the environment must be 
modified in the following way 

7<3-10~ 8 . (57) 

We also study the dependence of the decoherence time 



200 
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20 
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0.80 



T 



0.82 



0.84 



0.86 



0.88 



0.90 



FIG. 3. [Color Online] Dependence of the decoherence time 
as a function of the minimum velocity « m in- 

with the shape of the velocity profile. To do that we 
plot the expression in Eq. (56) as function of v min . 
As explained in Ref. [8] this magnitude must be near 
0.832-7r/T to perform the measurement. As observed in 
Fig. 3, the decoherence time varies smoothly with the 
velocity, therefore small deviations over 0.83 • 2tt/T do 
not present inconvenience with respect to decoherence as 
long as 7 < 10~ 8 and T < 100 • T H . 



B. Behavior of the Decoherence Time with 
Temperature. 

In this section we will study how temperature influ- 
ences the previous results. To achieve this we need to 
calculate the diffusion coefficients again, in the presence 
of nonzero temperature 



D(t) = 



I 2 



dv I dsv coth cos i^s cos ws. (58) 

Jo 2 



To calculate this integral in the low temperature limit we 
will use the following approximation 



100 



10 



0.1 



t_D 






T 






r~io -8 




y~10 -7 ■ 


T 









20 



40 



60 



80 



100 



FIG. 4. [Color Online] Dependence of the decoherence time 
for different relative noises in the force as a function of tem- 
perature. 



< 1 and B^L > i. The details of the calculation 
can be found in appendix A. The result for the integral 
of the diffusion coefficient to lowest order in temperature 
is given by 



, ,-./ on -qWtt _ 9 2(1 — costL^u) 
dtD{t,l3) = 1 2 —t D + 1 2 



, _ 9 W7T 

^7 2 -r*u 



2w 2 /3 2 /i 2 

27 2 
w 2 /3 2 ft 2 ' 



(60) 



Using this expansion we can conclude that the decoher- 
ence time at finite temperatures is 



t D (T ) 



2h 2 8k 2 B T 2 
r y 2 AvS 2 TTUjp 2 V uj 3 irh 2 



+ 0(7$). (61) 



Taking advantage of the expression in Eq. (61), lets 
also study the dependence of the decoherence time with 
the initial temperature of the field To The result is 
plotted in the Fig. 4. One can observe in the figure that 
the decoherence time does not change substantially with 
temperature for small enough coupling and therefore this 
makes no restriction on the temperature. Decoherence 
aside, taking into account the results obtained in Ref. 
[8], the entanglement that one would wish to measure 
as a signature for the Hawking effect is present only for 
To < 100 • Tff. It is also important to notice that the 
low temperature approximation is valid in the range T < 
100T H . 

Here we assumed that the system is initially in thermal 
equilibrium with respect to the environment. However, 
if the bath is at a higher temperature than the system 
one would expect that the bath heats up the system and 



coth 




si 



fitly 
1 

t3hv 
2 



< 1 
> 1. 



(59) 



and perform the integrals in the two regimes, both 



1 In general, To is the initial temperature of the environment. 
Since we consider that initially the system is in thermal equi- 
ibrium this is also the initial temperature of the system. 
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so washes out the Hawking signal. This featured is re- 
flected in the fact that the decoherence time decreases 
for increases temperatures T , which must now be inter- 
preted as the initial temperature of the bath. 



include all the relevant information regarding the envi- 
ronment. In the case of harmonic oscillators with a spec- 
tral density their expression was already given in Eqs. 
(24) and (25). 



V. DYNAMICS OF THE ENTANGLEMENT 
THROUGH CORRELATIONS 

In this section we will study the signature of the Hawk- 
ing effect through the shape of the correlation for a 
generic acoustic black hole, following the presentation 
given in Ref. [11]. 

Hawking radiation can be understood as a pair produc- 
tion of virtual particles, one of which falls into the black 
hole and the other, outgoing, becomes real, building up 
the Hawking radiation [1], as can be seen from Eq. (10). 
The role of the correlations between this pair of particles 
has been studied as a signature of the quantum Hawk- 
ing effect, see for example Refs. [8, 9, 13]. It was found 
that the entanglement between this pair is translated to 
a sharp peak of (Hl{xi, t)TlL{x, t)} as a function of x, 
where x\ is inside the black hole and x is outside, n is the 
canonical momentum conjugate to </>; the subscript L cor- 
responds to the left-moving modes. This magnitude was 
calculated using the Israel-Hartlc-Hawking state [13, 21- 
23] and it was also calculated numerically in the case of 
the circular ion trap [8] and BEC [9] . The result is shown 
in [13], Fig. 1. 

The formalism used in Ref. [13] does not seem to be 
useful in the case of an acoustic black hole as an open 
quantum system. The presence of the environment does 
not seem to be easily included. Moreover, in the com- 
putation a model of eternal black hole is used. There- 
fore, the inclusion of an initial thermal state and the time 
evolution of the correlation function cannot be obtained 
within previous procedures. 

In this section we want to resolve this issues to some 
extent. We will analyze the complete evolution starting 
from the flat background coupled to a massless scalar par- 
ticle, initially in a thermal state of temperature To, and 
"collapsing" toward an acoustic black hole background 
with an event horizon. 

To find the correlation functions in the quantum 
regime, we use the stochastic formalism developed by 
Calzetta et al in Ref. [24], reviewed in Ref. [19] and 
applied in Ref. [11] to this problem. This procedure will 
allow us to find the correlation functions through the 
semiclassical solutions of the Langevin EOM. We will as- 
sume that we study a generic acoustic black hole with 
a generic environment. The coarse grained effective ac- 
tion of a generic field, that emerges for example from the 
presence of an environment, can be cast into the form 

S cS [0, 4>'} = So W\ - So W] + Sif [4>, 4>'\ , (62) 



A. Stochastic Description 

The purely imaginary term of the effective action can 
be written in the following way 2 

J X>£ P[f]e _i *- = e***- N -.-'*-' , (63) 

where, for example <f> x = 4>{x) and repeated indexes are 
integrated. The probability distribution is gaussian 

P[£] ee N ( .e~^° ,n ^ A *'' , (64) 

where is a normalization constant, see Ref. [18, 19, 
24]. Therefore, the original action is the same as a 
stochastic process described by an action of the following 
type 

Sty, 0'] = Soft] - S [<//] + j-n^tfe - (65) 

We will start with general results regarding the correla- 
tion function and later use it to the particular case of 
acoustic black holes. 

To find the correlation functions we must solve the 
semiclassical EOMs for a realization of the stochastic 
force £(t, x) and then integrate over the possible solu- 
tions with a weight imposed by the probability distribu- 
tion P[£] and also integrate over the initial conditions. It 
is in the latter integration that the initial temperature of 
the field T makes its appearence. 

To start, we will make a distinction between the ho- 
mogeneous and inhomogeneous part of the solution, i.e. 

4>^(x, t) = <po(x,t)+ / dsdx'G Te t(t, s\X) x')£(s, x 1 ) (66) 
Jo 

where 4>o{x,t) is the noiseless solution satisfying 

a<p (t,x)+ [ dsD(t, s)<j) (s,x) = 0, (67) 
Jo 

where □ = V^V^ 1 and G rc t is the retarded Green function 
associated with the equation, satisfying 

□ G rc t+ / tft'D(t,t')Gret(*',s) =S(t-s)S(x'-x). (68) 
Jo 

We take the usual Painleve-Gullstrand-Lemaitre metric 
associated to a generic acoustic black hole 

ds 2 = -dt 2 + (dx - v(x, t)dt) 2 

= -(1 - v(x, t) 2 )dt 2 - 2v(x, t)dxdt + dx 2 , (69) 



where the free action is given by Eq. (12), the influence 
action by Eq. (23) and the kernels T)(x, x') and N(a;, x') 



2 In this section, to simplify the calculation we use c = h = 1. 
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where v(x, t) is the velocity profile, in principle arbitrary 
as long as it has a supersonic subsonic transition, i.e. an 
event horizon. In this section we now consider its time 
dependence instead. 

In this paper we work under the weak coupling approx- 
imation, in such a way that we can estimate 4>o beginning 
with the zeroth order limit of the EOM 0<fi ~ which 
explicitly, replacing the metric of Eq. (69), is given by 

[(d t + d x v(x,t))(d t +v(x,t)d x )-d*}<t>o(x,t)=0, (70) 

and moreover, we can approximate G re t from 

□G ret (*, s\x, x') ~ S(t - s)S(x' - x). (71) 

Therefore to first order in the coupling with the environ- 
ment we can consider the dissipation term in Eq. (67) as 
an inhomogeneity, replacing the complete solution <po by 
the solution of the Eq. (70), 



0<f>~£(x,t)- [ dt'T>(t,t')</>o(x,t'), 
Jo 



(72) 



which in turn has the solution 



<f>t(x,t) ~ 4>o(x,t) + / dsdx'G ret (t,s\x,x')£f(s,x'), 
Jo 

(73) 

where we define as 



dsT)(t, s)4>o{x, s) 



(74) 



Assuming that we have this solutions, the correlation 
function relevant to the study of the entanglement be- 
tween the Hawking pair phonons can be obtained from 



J/in' 



(75) 



where ( . . . ) in represents the mean value integrated over 
the initial conditions weighted with the thermal distribu- 
tion with To. To obtain the correlation function between 
the conjugate momentums one simply has to derive the 
equation above following the definition below 



IL(x,t) 



dej) d(f) 
dt+^dx- 



(76) 



Therefore, following the formalism presented in [24], 
<<M 2 > = ((0oi + G lx ^o2 + G 2x tl)) s . n 

= (4>Ol4>02)in ~ G2x>'D x i x ((f>Ox4'02)in 
— G 2x '~D x i x {4>Ox4>01 )in 

+GW<&'£ E > £ , in G 2 x. (77) 

Using the properties of the stochastic force probability 
distribution we get 

(0102)= (0O10O2)in - G 2x 'T) x ' x {(po X (p02}in 

— G 2x >T} x 'x(4>Ox4>Ol)in + Gi x >N x > x G2x- (78) 



In the case without environment, where the field <f> rep- 
resents particles in an acoustic black hole background, 
the previous formalism can be used but only the first 
term contributes. Therefore, in the general case one can 
write 



'192)0 



'\V2lc 



•0(A), 



(79) 



where A is the coupling constant with the environment. 
The subindices C and O means that it corresponds to the 
closed and open system, respectively. This simple case 
dropping the 0(A) terms already presents the main diffi- 
culties inherent with the calculation and we will develop 
a technique to deal with the equations in the following 
sections. 



B. Modes for the Wave Equation. 

To solve the EOM for the field we will use the method 
of characteristics, well known from mechanics of com- 
pressible fluids, see Ref. [25]. If one has a general differ- 
ential equation of the form 



Qu du 
a(x,t)- — h b(x,t) — + c(x,t)u 
ox ot 



0. 



(80) 



with initial condition u{x, 0) = f(x) then it can be solved 
in the following way. First find the characteristic curves, 
defined as 



dx 
ds 



a(x, t) 



and 



dt 
ds 



b(x,t). (81) 



Then find the evolution of u(x,t) along the characteris- 
tic curves finding the solution for the following ordinary 
differential equation 



du 
ds 



(x(s),t(s)) + c(x(s),t(s))u(x(s), t(s)) = 0. (82) 



Finally, having the congruence of characteristic curves, to 
know u(x, t) invert (x, t) <-> (xq, s) where Xq is the initial 
condition and s the parameter along the characteristic 
that goes through (x,t). Then, u(x,t) = /(xq). 

In the rest of this section we will use this method to 
find the solution of the EOM for the field. 

Before attempting to solve for the modes we have to 
give a specific velocity profile, 

Vmin para — oo < x < —a 
v(x,t) = u{t) x ^ 1 + kx para — a < x < a (83) 
fmax para a < x < oo, 

where a(t) is a function that must satisfy er(0) = 0, to 
guaranty that at t = the metric is flat and a(t 
t) = 1, in order for the background to "collapse" in an 
stable acoustic black hole after the time interval t. A 
particular function that satisfies this requirement and is 
easy to manipulate is 



a(t) = tanh 



t 



(84) 
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We will solve the equations first for the region —a < 
x < a without taking into account the rest of the space. 
These solutions are the same in the other regions upon 
the following changes k = and a i— > v max /mm <T - After 
this, we will see how to introduce them properly and thus 
obtain the full solution. 

Nevertheless, the problem that the presence of differ- 
ent regions introduce are the interfaces. For example, a 
characteristic curve that starts in —a < xq < a eventu- 
ally reach the point x(si) = a and for s > Sj the char- 
acteristic to use is the one corresponding to the region 
a < x < oo. In this section we will see that the solutions 
that do not go through this boundaries do not present the 
characteristic peak associated with the entanglement of 
the Hawking pair. In the following section we will study 
the full solution and we will learn how the entanglement 
is developed and its relationship to this issues. 

The equation we have to solve is Eq. (70) which can 
be cast into the following form 

(d t +d x v(x,t)+d x )(d t +v(x,t)d x -d x )<p (x,t) = 0. (85) 



Defining the operators 

dL = dt 
d R = 



v(x, t) 



dx 



d_ 
dt 



d_ 
dx 

d , , d 
dx VM+ dx- 



(86) 
(87) 



the equation can be written as d^d^ — 0. One can first 
solve <9l</>l = 0. Then one has to solve Br^r = and 
finally 8l4>r = <pR m order to obtain (f>R. This way, the 
more general solution is 00 — 4>l + 4>R- Both left and 
right moving components can not be solved separately 
since [d L ,d R ] ^ 0. 

The procedure to find the solution with the character- 
istic curves is developed in appendix B. For example, the 
characteristic curves for the left moving modes are given 
by 



x(t) =exp k / a(s)ds 



x x 



To make contact with the usual expansion in flat field 
theory we cast the solution into the following form 



x) = dk(u k (x)a k + H.c), 



(89) 



where H.c. means taking the hcrmitian conjugate of the 
expression and the modes are given by 



1 - e(fc)2i|fc| / ds 
Jo 



The coefficients a k and a\ corresponds at early times 
to the creation and annihilation operators upon quanti- 
zation and they would carry the subindex "in" . 

To compute the correlation function when an environ- 
ment is added we need the retarded Green function that 
we compute as G iet (t,t'\x, x') — [4>(x, t), <j>(x', t')]Q(t — 
f). 



C. Entanglement: Closed System. 

In this section we study the behavior of the two-point 
function of the left moving part of the momentum, n^. 
Using the expansion given in Eq. (89) and the definition 
of the left moving part of the momentum, Eq. (76), we 
can obtain the correlation function after integrating over 
the initial conditions 



(n L { Xl ,t)ii L (x 2 ,t)) 



dk 



2k 



-ifc(xi-£c 2 )e~" J o J* a (s)ds 



coth 



(91) 



with the usual definition p = (/c^To) -1 . 

Either for the region —a < x < a where the velocity 
changes with position, or the regions where the velocity 
is homogeneous (as stated, upon the replacements k h > 
and a t-¥ va), the expression below depends only on 
(xi — X2). Therefore, this expression cannot present the 
characteristic peak discussed. Accordingly, the regions 
where this solution is valid do not present the signature 
of the quantum Hawking effect, as explained in previous 
sections. 




FIG. 5. Schematic illustration of a characteristic curve that 
intersects the interface in x = a and starts at xo. 

The correct way of getting the solutions of different 
regions is matching the characteristic curves as shown 
in Fig. 5. For example, if x > a and different regions 
share the characteristic then xq is obtained matching the 
characteristic for v{x,t) constant with the characteristic 
for x < a and finally finding the Xq corresponding to 
the latter. The signature of the entanglement between 
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the Hawking pair was shown to be present in the left 
component of the correlation so we are only interested in 
the curves associated with this modes. 

To carry on this procedure, we approximate for long 
times (t t) the characteristic for |x| < a as x(t) ~ 
Xoe Kt . In the region x > a, the characteristic is 

x(t) =A-t+ /(t)tw, (92) 

where A is a constant of integration, f(t) = J Q dsa(s) 
and v m ax ~ 1 + Kffl. The relationship between A and xo 
can be found recalling that 

x(t a ) = a =4- t a = k" 1 loga/x 0l (93) 

because the |ac| < a characteristic reach x — a at t a and 

A - t a + f(t a )v max = a, (94) 

because the x > a characteristic must match the previous 
one at x(t a ) = a- 

From this two equations, for times t > t, we can ob- 
tain A as a function of xq and then use xq as a function 
of x, t and put it in the e tkx ° factor, finding 



xo = ae^ x+t -f^ 



- a ^ a , for x > a and (95) 
x = - ae - (x+t - f{t)v ^- a)/a 1 for x < -a. (96) 

Therefore the left moving modes are given by 
1 




, for x > a 



e lKae , for x < -a. 



Finally, we calculate (<fii,(xi,t)(j>L(x2,t)) for an initial 



1 <n L (-v) n t (jt)> | 
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FIG. 6. Plot of the absolute value of the left-moving compo- 
nent of the correlation function for long times t — lOOr. The 
maximum velocity is 1.1, the minimum 0.9 and the initial 
temperature is zero To = 0. 

thermal state at temperature /3 = (fc^To) -1 , and deriv- 
ing with respect to x,t according to Eq. (76) in order 



to obtain the momentum two point function. The final 
result is 

(U L ( Xl ,t)U L (x 2 ,t)) = -^e 5 ^[cosh(^)]"^ X 

cosechr^-e ^^coshf-)! " 
a 1 v t /j 



x (e~ 



e 2 [cosh(^)]^)]^ , (97) 



where A/Ew = w ma x — / + v m i n . This dependence of the 
correlation function between y = —x% — 4, 6 and x 2 = x 
is shown in Fig. 6. If y is far from the event horizon, 
placed at x = 0, the peak in the correlation will be far 
of the event horizon too, in the inverse direction. Even 
though the peak does not seems to modify with increasing 
y, it disappear for y big enough, as we will see below. 
This can be concluded directly from the characteristic 
curves. If one follows the characteristic with xq = a then 
everything to the right will not be aware that there is an 
event horizon in the space; therefore, there should not 
be entanglement (peak in correlation) and there is none 
indeed. This is because in this region we have to use 
the solutions given by Eq. (91), which was shown not to 
present the signature. 

Moreover, we can compare Fig. 6 of this paper with 
Fig. 1 of reference [13]. We can conclude that for long 
times f > t, the one calculated for the black hole with 
the Israel-Hartle-Hawking coincide with the calculation 
in the present paper and in Ref. [11], that includes the 
collapse. Their behavior is the same, but little differences 
may arise from the fact that the velocity profile as a func- 
tion of x is different, since they use a smooth tanh(a:/a) 
profile. 

Nonetheless, this formalism allows us also to study the 
evolution of the correlation function. As we explained, 
there are two special characteristic curves, that start 
at xo — ±a, and they separate two regions, one that 
presents the peak that reflects the entanglement and an- 
other one that does not. Therefore, the interface that 
separates both regions, which is shown in Fig. 7, is given 
by 



x±(t) = ±a(l + kJ a(t)di 



= ±a 1 + kt In cosh 



(98) 



As can be seen in Fig. 7(b) if we study the correlation at 
fixed position x there exist a definite time tE such that 
for t < tE there is no entanglement and it is generated 
later for t > tE- This time can be found solving the 
equation x+(t£;) = x. 

In turn, if we study the correlation at fixed time t, then 
there exists a position xe such that if a < x < xe there 
is no peak in the correlation and for x > xe it is present. 

For example, for fixed time t = lOOr with the param- 
eters of Fig. 6 one has the presence of the correlation 
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FIG. 7. [Color Online] (a) Illustration of the regions where the 
peak of the correlation function is present, and therefore the 
entanglement between the Hawking pair particles and those 
that not. (b) Behavior for fixed time or position. Red: no 
correlation peak. Blue: presence of correlation peak. 



peak only for \x\ < 11, therefore it is correct to use the 
expression in Eq. (97) and not that of Eq. (91) to make 
the plot in Fig. 6. The fact that this transitions are not 
smooth is due to the abrupt changes at x = ±a in the 
velocity profile. 

These two features of the solution found can be com- 
pare with the numerical simulations done for the ion ring, 
[8] . The fact that the correlation peak has a finite spatial 
extension can be clearly seen in Fig. 9 of Ref . [8] , where 
this magnitude is calculated numerically. The fact the 
the entanglement needs a definite time to develop is also 
observed in Fig. 14 of Ref. [8] through the numerical 
computation of the logarithmic negativity. 

This formalism also allows us to study how the corre- 
lation, and therefore the entanglement, depends on the 
initial temperature of the field T , as can be seen in Eq. 
(97). In particular, our results should reproduce the be- 
havior shown in Fig. 14 of Ref. [8] , and the loss of entan- 
glement of the Hawking pair at high temperatures. First 
we compute the Hawking temperature, given by (see Ref. 
[12]) 



T, 



H 



i k 



2tt 



2a 



(99) 



In Fig. 8 we plot the dependence of the peak in the 
correlation as a function of the temperature in units of 
the Hawking temperature. We observe how the peak is 



0.25- 

(n L (-l()a, 10-r)n,]x, ltf^T 
20 Black Hole 

0.15 




FIG. 8. Plot of the absolute value of the correlation function 
of the momentum for long times t — lOOr. The maximum ve- 
locity is 1.1, the minimum is 0.9 and the initial temperatures 
are T = 0, 20 • T H ,60 • T H . 



diluted and therefore how the Hawking pair entanglement 
is lost for high enough initial temperatures. 



D. Entanglement: Open System. 

In this last section, we will try to estimate the correc- 
tion introduced by the presence of an environment in the 
correlation function. Therefore, we will use the expres- 
sion presented previously 



n<P2 



- (<t>Ol4>02) in — G2 x '"D x i x {(f)Ox<t > 02)m 
-G2x''Dx'x(4 l Ox4>Ol)in + G\ x ' N^j G2X , (100) 



together with the solution for the modes found in the 
previous sections. Regarding the environment, we will 
use an ohmic bosonic bath with bilinear coupling. The 
dissipation kernel is given by 



D( s ,,s') = -A 2 |-<5( 
as 



(101) 



In this case, obtaining an expression analytic for 
(HlIIl) is more involved since the integrals must be per- 
formed over all possible regions and the final dependence 
of k make the integrals more complicated. To solve the 
second inconvenience we will study each mode separately, 
characterized by k. We will use the expression in Eq. 
(101) to eliminate the integrals in s' , after integrating by 
parts the Dirac delta function. The remaining integrals 
can be approximated by the main contribution in the 
coincidence limit. Regarding the last term, for high tem- 
peratures the noise kernel is N(si,S2) = X 2 TqS(si — S2) 
and the resulting integral is done using basic properties 
of Green functions. After performing the derivatives nec- 
essary to obtain (H^Hl) we can estimate the correlation 
function at high temperatures. Fig. 9 shows the result of 
the relative contribution of the term associated with the 
environment with respect to the solution for the closed 
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FIG. 9. Relative contribution e r for high temperatures, 
T ~ 100 ■ Th, evaluated on the position of the peak in the 
correlation function and for coupling A ~ 10~ 7 . 



black hole, given by 

(n L1 n £2 ) c (fc) - (n L1 n L2 ) (fc) 



<n £1 n i2 ) c (fc) 

where the k dependence means that 



(n L1 n 



L2 



dk(n L1 u L2 )(k). 



(102) 



(103) 



This magnitude was computed for small k since it de- 
creases with increasing k. For small times the influence 
of the environment is negligible and the system behaves 
as an effectively closed system. For later times the rela- 
tive contribution begins to grow. Even though we could 
not integrate to obtain the result as a function of x, t 
we expect that this abrupt increase of the environmen- 
tal influence produces decoherence and we take it as a 
signature of the loss of correlation and therefore, the en- 
tanglement between the Hawking pair particles. 



VI. CLOSING REMARKS 

In this paper we studied the influence of an environ- 
ment in acoustic black holes as an open quantum sys- 
tem. We used the circular ion trap setup but the anal- 
ysis can be applied to any acoustic black hole, only de- 
tails relevant to the environment could change. We put 
our attention to the process of decoherence induced by 
this environment in contradistinction with the quantum 
to classical transition considered previously induced by 
thermal equilibrium effects [8]. As we explained in the 
Introduction, since the Hawking effect is purely quan- 
tum mechanical, then a quantum to classical transition 
would jeopardize the possible measurement of the radi- 
ation, therefore decoherence must be controlled in the 
experiment. 

We used the open quantum system approach to cal- 
culate de decoherence time of an acoustic black hole in 



presence of an environment for both zero and finite (but 
small) temperatures [11]. Taking into account that the 
decoherence time is smaller than the approximate mea- 
surement time and worst, the "collapse" time, improved 
parameters are needed and provided here to make the 
measurement possible. This analysis could be easily ex- 
tended to other realizations of acoustic black hole such 
as BEC, moving dielectrics, waveguides, etc. 

In order to achieve generality for this work to be in- 
teresting for any acoustic black hole, we made certain 
approximations. It would be interesting for future works 
to perform this derivation numerically taking as a start- 
ing point the exact dynamics of the particular acoustic 
black holes' setup. Moreover, regarding the environment 
one could consider more involved bosonic reservoirs that 
for example gets dragged to some extent by the moving 
ions instead of being fixed with respect to the laboratory. 

We also provide a derivation of the correlation between 
the Hawking phonons as a function of time and temper- 
ature and we check that the relative contribution of the 
environment to this magnitude is irrelevant. The deriva- 
tion included the transient due to the "collapse" of the 
acoustic black hole and could be applied to the usual 
Schwarchild black hole as in [22]. 

The main complication in this approach is the presence 
of boundaries between different velocity profiles, since 
the overall is not continuously diffcrcntiable. For future 
works, it would be interesting to find a smooth spatial 
velocity profile that presents an event horizon and allows 
to carry an analytic treatment to the end 3 , in order to 
obtain the correlation in the presence of an environment 
as a function of t and x instead of Fig. 9. 

Finally, another interesting extension of this work 
would be to study the quantum to classical transition due 
to the non linearities in the original hamiltonian, Eq. (1). 
The non equilibrium self interaction of the degrees of free- 
dom in this case (as opposed to the thermal equilibrium 
situation considered in Ref. [8]) also introduces dissipa- 
tion and noise in the same fashion as an environment and 
also induces decoherence. The decoherence time associ- 
ated with this effect depends on intrinsic parameters of 
the system such as the mass, electric charge of the ions, 
etc.; as opposed to the external parameter 7 considered 
here. Therefore, this study would impose critical bounds 
on the realizability of the measurement of the Hawking 
radiation with any setup. 



ACKNOWLEDGEMENTS 

This work was supported by UBA, CONICET and AN- 
PCyT. GJT is also supported by CIN. The authors thank 
Esteban Calzetta, Juan Ignacio Cirac, Diana Lopez Nacir 
and Juan Pablo Paz for useful discussions. 



3 For example, we tried to use a tanhx/a profile but exact analyt- 
ical expressions for the characteristic curves could not be found. 
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Appendix A: Diffusion coefficient at finite 
temperature 

The dissipation kernel for finite temperature is denned 



as 



7 2 f°° f* j3hv 
D(t) = — / dv \ dsv coth — — cos vs cos los. (Al) 
2 Jo Jo 2 

As stated in section IV B, we will use the following ap- 
proximation 



coth 



f3Hv 



/3hv bl 2 ^ 1 

si Bg- » 1. 



1 



(A2) 



Having this in mind, the coefficient can be approximated 
integrating separately over the regions where the different 
approximations are valid, i.e. 



dv 



(i 



209ft)" 1 



dv v coth 



x / ds cos los cos vs 
Jo 

2(f)K)- 1 2 ,t 

dv— I ds cos los cos vs 
fib Jo 

dvv I ds cos lus cos vs. 

2(fih)- 1 Jo 



The first integral is 

2(/3h)- 1 2 ,t 

dv— / ds cos lus cos vs = 
fib Jo 

For small 2(j3h)~ 1 we can approximate 

4 sin Lot 



(A3) 



4t 



and therefore 



"^TT^ / us cos ws cos j/s — . (A4) 
P" Jo 



Regarding the second term, we will use the following 



dvv / ds cos uis cos vs = 

2(fih)- 1 Jo 

f* 

dvv I ds cos uis cos vs 

n Jo 

2(pny 1 rt 

dvv I ds cos los cos vs. (A5) 

n Jo 



The first term was calculated for the Tq = case and 
with the result LoSi(cot). The second term is 



2 sin Lot 



203ft)- 1 rt 

dvv / cis cos los cos vs « — . (A6) 



Finally the new diffusion coefficient is given by 

2T 2 D(t, P) = LoSKLot) + 2 -^§ + 0(r 3 )- (A7) 

When one takes lo = and t = the correction to order 
~ 0(Tq) vanishes. Therefore, a thermal state of finite 
temperature does not affect the relationship between the 
parameters 7 y 7. 

The next step is just to integrate the Eq. (A7) in order 
to obtain the relevant diffusion coefficient. 



Appendix B: Modes of the Wave Equation 

The equation we have to solve is Eq. (70) which as 
noted above, can be written as 9r<9l0 = 0. This way, 
the more general solution is <j>o = 4>l + <f>R, in terms of 
the left and right moving modes. 

The first equation to solve is <9l</>l =0 



_ + („(*,*) -1)_ = 0. 



(Bl) 



This equation is of the same form that the Eq. (80) with 
b(x, t) = 1, a(x, t) = v(x, i) — 1 = a{t)nx + <r(t) — 1 and 
c(x,t) = 0. 

The definition of the characteristic curves is 



dx 
ds 



a{t)nx + <j(t) — 1 and 



dt 
ds 



= 1. (B2) 



From the second equation one obtains t = s, where we 
impose that t(0) = 0. Using this relationship in the first 
equation we get 



dx 
~dt 



a(t)nx — a(t) — 1. 



(B3) 



The solution of this equation is 

x(t) = e K f<X s ) ds 

x - [ {l-aitye-Kfo^'^s'ds), (B4) 



where we used the condition x(0) = x . Then, t(xo, s) = 
t = s and inverting to get xq{x, i) we obtain 

x {x,t) =xe- K ti a ^ ds 

+ f (l-a(t))e- K J'°' 7 ( s ') ds 'ds. (B5) 



The differential equation to solve for the evolution of the 
field along the curve is simply 



ds 



=> <f>(x(xo,s),t(xo,s)) = 0(xo,O). (B6) 
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Finally the solution is 

+ (l-a(t))e- K J'o^ s '^ s 'd S ,oj. (B7) 

To rewrite it in a useful way, taking into account that 
we will have to eventually integrate over the initial con- 
ditions, we expand in modes <f>(x, 0) = J dk/2Tre lkx <j)k, 
then the solution for the field is 

H(x,t) = J ^ fc exp{ifc(se-"/ M«)* 

+ J\l - a{t))e~ K fo ^ s ') ds 'ds^ | . (B8) 

Having an expression for the left moving modes, we 
now have to solve the equation for the right modes, 
d R (j)R = 0, 



d d , , 



+ a(t)n + cr(£)(l + (tx)— + 



-) 




dxj 


4>R 


~) 




dxj 


4>R 


~) 




dxj 


4>R 



IS 



d 
dt 

d_ 
Ft 



Again, t = s, but the spatial part of the characteristic 

x(t) = e K K a(s)ds 

(x + j\o{s) + l)e" K /o °(>')*>'d8^ , (BIO) 

and after taking the inverse xo 

x (x,t)=xe- K ti^ ds 

- f (l + a(s))e- K fo^ s ') ds 'ds. (Bll) 
Jo 

Now the equation for the value <f)R along the character- 
istics is 



dt 



= —a(t)K<j>, 



(B12) 



then 



<l>(x(x ,8),t(x ,8)) = cfr(xo,0)e- K f°^ ds '. (B13) 
Therefore, expanding again in modes 

- jf (l + (r(t))e- ,, J'o'(«')*>'d^ 

-k^ cr(s)ds|. (B14) 



Taking advantage of this result for <pR, the right moving 
mode can be found from (dt + v(x,t)d x — d x )<j)R — 4>r. 
To solve this equation we propose a solution of the form 



<t> R {x,t)= y ^ k {x,t)e^ik(xe- K ti°^ ds 

+ y (l- ( r( S ))e- K ^ CT(s ' )ds 'rf^|. (B15) 



If we use this result in the equation and use that the 
exponential is a solution of the homogeneous equation of 
the left moving modes, we obtain 



>L<?R 



(Ik 
■t 



d L (f>k(x,t)exp 



ik(xe~ K ^ (y ^ ds 



+ 



f {l-a{s))e- K fo^') ds 'ds) 
Jo 



(B16) 



On the other hand, the solution <f>R can be rewritten as 
4> R (x,t)= J ^-4> k exp[-2ik J e - K f° u ^ ds 'ds- k J a(s)ds} 

xe ifc(xe-J3"C> , "+/ '(l-aW)e-»J"tf''<''""'«fa) > ( B1? ) 

Equating both sides of this equation gives 

d L <t> k (x,t) = ~ 4>ke -^S t a e-!S^'^'ds-.S^{s)ds_ (Blg) 

To find the particular solution we can use a 4>k such that 
4>k(x,t) = &k(t) and therefore 

dt<l> k (t) = fa-™ ft *>-"ti "M* (B19) 

Integrating in time, the solution is 

ik(xe- K So "( s ' ds +/ t (l-<T( S ))e- K Jo <*(«')<«•' ds ) 

/ 2tt fc 

f dse- 2ik /o /o "' " < '" ) '"" /o . (B20) 

Jo 



t(*>*)= / 2 _ <l '- 



The remaining steps are just to write down the full so- 
lution and cast it into the usual mode expansion form. 
To do this one has to write and 4>kR hi terms of the 
Fourier transform of the initial condition of the field and 
the conjugate momentum <fik{t = 0) and Hk(t = 0). 
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